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• Premise of the study: Microsatellite loci were developed for the rare, Himalayan, endemic haploid lichen fungus, Lobaria 
pindarensis, to study its population subdivision and the species' response to forest disturbance and fragmentation. 

• Methods and Results: We developed 1 8 polymorphic microsatellite markers using 454 pyrosequencing data and assessed them 
in 109 individuals. The number of alleles per locus ranged from three to 1 1 with an average of 6.9. Nei's unbiased gene diver- 
sity, averaged over loci, ranged from 0.514 to 0.685 in the three populations studied. The cross-amplification success with re- 
lated species (L. chinensis, L. gyrophorica, L. isidiophora, L. orientalis, L. pulmonaria, L. spathulata, and Lobaria sp.) was 
generally high and decreased with decreasing relationship to L. pindarensis. 

• Conclusions: The new markers will allow the study of genetic diversity and differentiation within L. pindarensis across its 
distribution. Moreover, they will enable us to study the effects of forest management on the genetic population structure of this 
tree-colonizing lichen and to carry out population genetic studies of related species in East Asia. 
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Lobaria pindarensis Rasanen (Lobariaceae, Peltigerales) is a 
foliose lichen species known from mountain forests and open 
woodlands in the Himalayas of Bhutan, India, and Nepal. The 
lichen thallus is haploid and it mainly disperses with vegetative 
propagules, but sexual reproduction with ascospores can also 
occur (Scheidegger et al., 2010). The lichen disperses locally, 
thus sharing ecological traits with L. pulmonaria (L.) Hoffm. 
(Scheidegger and Werth, 2009; Scheidegger et al., 2012). Al- 
though microsatellite markers are available for L. pulmonaria 
(Dal Grande et al., 2010; Werth et al., 2013), only three markers 
(LPu32425, LPu40211, and LPu34888) published by Werth 
et al. (2013) reveal small, multiple bands when amplified with 
L. pindarensis. All other published markers do not amplify with 
L. pindarensis. Here, we develop microsatellite markers to 
study the impact of land use and habitat fragmentation on gene 
flow of this dispersal-limited lichen (Scheidegger et al., 2010). 
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METHODS AND RESULTS 

Ten specimens of L. pindarensis, collected in two valleys in Nepal (Table 1 ; 
Manaslu Conservation Area [MCA] and Sagarmatha National Park [SNP]), 
were used for total DNA extraction with the QIAGEN Plant Mini Kit 
(QIAGEN, Hilden, Germany). Subsequently, whole genome 454 pyrosequencing 
of pooled DNA samples was performed using a Roche GS FLX sequencer 
to generate a sufficient number of microsatellite loci. Library preparation 
and sequencing were performed by Microsynth (Balgach, Switzerland). 
Shotgun libraries were prepared using the GS FLX Titanium Rapid Library 
Preparation Kit (Roche Diagnostics, Basel, Switzerland), while Microsynth 
provided barcode adapters. Out of a l/4th run, we obtained 233,260 reads 
of an average length of 314 bases for a total of 73,191,881 bases. The unas- 
sembled sequences were screened for all possible sequence motifs of di-, 
tri-, tetra-, and pentanucleotide microsatellites using Primer3, as imple- 
mented in MSATCOMMANDER version 1.0.2 alpha (Rozen and Skaletsky, 
2000; Faircloth, 2008). Microsatellites with motifs repeated at least eight 
times (for dinucleotides) or six (for all others) were chosen. For each lo- 
cus, primer pairs were developed with MSATCOMMANDER using the de- 
fault parameters. 

Using all reads, MSATCOMMANDER found 1021 primer pairs that ful- 
filled the default primer parameters. Subsequently, 656 pairs were discarded 
either because they contained unfavorable secondary structure, primer dimer 
formation, or mononucleotide repeats in the flanking region, or because they 
were duplicates, which were detected after alignment using CLC DNA Work- 
bench 5 (CLC bio, Aarhus, Denmark). The remaining 365 sequences were 
verified one by one using ntBLAST with the megablast option (http://www. 
ncbi.nlm.nih.gov/blast) to exclude those that were highly similar to algae, 
plants, or microorganisms that are often present in environmental samples. To 
test for cross-amplification with the photobiont of L. pindarensis, Dictyochlo- 
ropsis reticulata (Tschermak-Woess) Tschermak-Woess, PCRs of the re- 
maining 116 primer pairs (including 44 di-, 65 tri-, and 7 tetranucleotides) 
were run using DNA from an axenic culture of D. reticulata (Dal Grande 
et al., 2010, 2012; Widmer et al., 2010, 2012). The 56 loci that produced posi- 
tive PCR reactions were excluded from further analyses because they were 
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Table 1. Characteristics of 18 polymorphic microsatellite loci developed 
for Lobaria pindarensis* and screened in 109 individuals. 



Locus 




Total 


MCA (n = 36) 


SNP (n = 43) 


KCA( 


^ = 30) 


N 


A 


A 


if 


A 


H t 


A 


H e 


LpiOl 


106 


4 


3 


0.643 


4 


0.615 


3 


0.587 


Lpi02 


109 


5 


4 


0.652 


4 


0.568 


4 


0.524 


Lpi03 


109 


5 


3 


0.160 


4 


0.295 


4 


0.582 


Lpi04 


108 


8 


6 


0.635 


7 


0.762 


6 


0.800 


Lpi05 


109 


7 


4 


0.162 


5 


0.666 


6 


0.715 


Lpi06 


109 


7 


4 


0.463 


5 


0.636 


6 


0.747 


Lpi07 


105 


9 


7 


0.567 




0.681 


5 


0.690 


Lpi08 


108 


5 


3 


0.565 


5 


0.741 


5 


0.594 


Lpi09 


109 


10 


8 


0.700 


5 


0.260 


5 


0.556 


LpilO 


109 


11 


5 


0.754 


8 


0.856 


8 


0.779 


Lpill 


108 


8 


5 


0.459 


3 


0.671 


6 


0.820 


Lpil2 


109 


4 


3 


0.256 


3 


0.456 


4 


0.724 


Lpil3 


109 


7 


7 


0.752 


5 


0.617 


7 


0.726 


Lpil4 


109 


10 


6 


0.308 


7 


0.780 


7 


0.786 


Lpil5 


109 


6 


4 


0.760 


5 


0.767 


5 


0.501 


Lpil6 


108 


10 


6 


0.816 


9 


0.791 


7 


0.788 


Lpil7 


96 


6 


4 


0.437 


5 


0.692 


6 


0.869 


Lpil8 


81 


3 


2 


0.170 


3 


0.462 


2 


0.533 


Average 




6.944 


4.667 0.514 


5.167 


0.629 


5.333 


0.685 



Note: A = number of alleles; H Q = Nei's unbiased gene diversity; n = 
number of samples per population; N = total number of samples analyzed. 

a Populations used in the study: MCA = Manaslu Conservation Area, Gorkha 
District (28°27.641'N, 85°2.803'E); SNP = Sagarmatha National Park, 
Solukhumbu District (27°48.871 / N, 86°43.016 , E); KCA = Kanchenjunga 
Conservation Area, Taplejung District (27°41.546 , N, 87°45.607 , E). Voucher 
specimens (collector numbers SD1164_IZ70ZO_131338/1_MCA, SD135_ 
IZ70ZO_131338/1_SNP, and SD268_IZ70ZO_131338/1_KCA) were 
collected in 201 1 and 2012 and were deposited at Tribhuvan University, 
Kirtipur, Nepal (TUCH), and in the frozen herbarium at the Swiss Federal 
Research Institute WSL, Birmensdorf, Switzerland. 



single products, all at an annealing temperature of 57°C. Polymorphism of the 
48 microsatellite loci was initially tested on a subset of seven individuals 
from three valleys, including three individuals from MCA, three from SNP, 
and two from Kanchenjunga Conservation Area (KCA), resulting in the de- 
tection of 25 polymorphic loci. These 25 loci were then tested in one popu- 
lation of 48 specimens from MCA, resulting in 18 loci with satisfactory 
amplification. 

To characterize these 18 polymorphic L. pindarensis loci, PCRs of 109 indi- 
viduals from three valleys were conducted (Table 1). The PCR protocol used 
fluorescent forward primers and the reaction was adjusted to the following con- 
ditions: 5 min at 95°C; followed by 25 cycles of 30 s at 95°C, 90 s at 57°C, and 
30 s at 72°C; with a final extension of 60 min at 60°C. All PCR products ob- 
tained were multiplexed (Table 2) and run on a 3130x1 DNA Analyzer with a 
GeneScan 500 LIZ Size Standard (G5 dye set) for fragment analysis (both by 
Life Technologies). Alleles were sized using GeneMapper version 3.7 (Life 
Technologies), and the variability of each microsatellite locus was measured by 
counting the number of alleles and calculating gene diversity using Arlequin 
version 3.11 (Excoffier et al., 2005). Trinucleotide microsatellites (n = 15) were 
the most common loci detected among the 18 microsatellite motifs (Table 2). 
The microsatellite loci produced 3-11 alleles per locus with an average of 6.9, 
and mean gene diversities over three populations varied from 0.514 to 0.685 
(Table 1). 

Cross-species amplification of seven closely related species of Lobaria 
(Schreb.) Hoffm. was tested on one specimen of each species (Appendix 1), 
applying the same PCR conditions established for L. pindarensis. All fragments 
were sequenced according to Cornejo and Scheidegger (2010) except LpiOl 
and Lpi05, which were verified on an agarose gel but not sequenced. The trans- 
ferability was high. Only one locus (Lpi05) did not amplify in L. chinensis 
Yoshim. However, several loci contained insertions within the flanking regions, 
and in others the microsatellite was disrupted (imperfect or interrupted mi- 
crosatellite). In some loci, the microsatellite sequences were reduced or dis- 
appeared completely, as in LpilO in L. isidiophora Yoshim. and Lpil6 in 
L. gyrophorica Yoshim., L. pulmonaria, L. spathulata (Inum.) Yoshim., 
and Lobaria sp. (Fig. 1). In general, the cross-amplification success of Lpi 
markers decreased with decreasing relationship to L. pindarensis, being lowest 
in L. pulmonaria and Lobaria sp. (with four and six loci not amplifying, 
respectively). 



considered alga-specific rather than fungus-specific. For PCR amplifica- 
tion, forward primers were labeled with an Ml 3 tag (5'-TGTAAAACGA- 
CGGCCAGT-3') (Schuelke, 2000). PCR reactions were performed in a total 
volume of 10 uL containing 1 uL of -1-5 ng genomic DNA, 0.5 uL of 5 uM 
forward and reverse primers, and lx Type-it Multiplex PCR Master Mix 
(QIAGEN). All PCRs were performed on Veriti Thermal Cyclers (Life Tech- 
nologies, Carlsbad, California, USA). The PCR reactions were assessed in a 
temperature gradient increasing by one-degree increments from 56-6 1°C, 
performed with the Type-it Microsatellite PCR Kit (QIAGEN) according to 
the manufacturer' s protocol, and under the following conditions: denaturation 
for 5 min at 95°C; followed by 33 cycles of 30 s at 95°C, 90 s at 56-6 1°C, and 
30 s at 72°C; then for the M13-tag binding an additional eight cycles of 30 s 
at 95°C, 90 s at 53°C, and 30 s at 72°C; with a final extension of 30 min 
at 60°C. 

The 60 remaining loci were assessed for amplification for the fungal 
component of L. pindarensis under the same conditions as above and using 
the total DNA of L. pindarensis. Out of these 60 loci, 48 produced specific 



CONCLUSIONS 

Fungi, algae, and/or cyanobacteria live in close contact 
within the lichen thallus and hence the manual separation of 
symbionts for later molecular analyses is technically unfeasi- 
ble. Therefore, symbiont- specific genetic markers have to be 
used in population genetic studies of lichens (Widmer et al., 
2010). The newly developed, highly variable fungus-specific 
markers reported here will allow detailed studies on regional 
genetic differentiation, effects of forest disturbance on ge- 
netic diversity, and the contributions of clonal and sexual 
reproduction in this lichen species. Moreover, the flanking re- 
gions of the microsatellites will be used for sequence analy- 
ses in future phylogenetic studies of related taxa of the genus 
Lobaria. 
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Fig. 1. Alignment of the Lpil6 sequence containing a trinucleotide microsatellite region. The flanking region was excluded from the graphics. This 
locus was initially developed for Lobaria pindarensis. The first four species contain a microsatellite with n > 9 repeats. The following two species have 
n = 3 repeats and are not considered microsatellites. Finally, in L. pulmonaria and Lobaria sp. this locus did not evolve a microsatellite sequence. 
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Table 2. Overview of the microsatellite loci developed for the lichen fungus Lobaria pindarensis. 



Allele size GenBank 



Locus 




Primer sequences (5 / -3 / ) 


T a (°C) 


Repeat motif 


Fluorescent dye 


range (bp) 


Multiplex 21 


accession no. 


LpiOl 


F: 
R: 


TTTGCGGTATAATCGACGCG 
CACACGACGTCACTGTCTTG 


57 
57 


(CGT) 10 


FAM 


255-264 


3 


KF318149 


Lpi02 


F: 
R: 


GGGATTGCAGAGGGATTTGC 
CATTTCCACTCCGTCAACCC 


57 
57 


(GADio 


VIC 


164-182 


2 


KF318150 


Lpi03 


F: 
R: 


CCCATTATGCCATGTCCTGC 
AGGAGGATAGTGATGGTCGC 


57 
57 


(CTT) 9 


FAM 


346-358 


2 


KF318151 


Lpi04 


F: 
R: 


CAGAACTAGGCCGCATTGTG 
TGAGCACGTTGTTTCACTCG 


57 
57 


(GTT) 10 


VIC 


89-122 


1 


KF318152 


Lpi05 


F: 
R: 


GACTGGCCGGCAATTAGTAG 
TGAAAGGGTCTTGTTGCGAAC 


57 
57 


(CTT) 9 


VIC 


111-154 


2 


KF318153 


Lpi06 


F: 
R: 


GCGTATTGGAGATGGCGATG 
GGCCTGAACATGGAATGCAG 


57 
57 


(GAT) 9 


PET 


148-194 


2 


KF318154 


Lpi07 


F: 
R: 


CAAGCCACCCACTCATTTCG 
GCTACACGTTTGGGTCTTCC 


57 
57 


(CTT) 10 


NED 


250-277 


1 


KF318155 


Lpi08 


F: 
R: 


CTTTCCTCGCTGCAGAACTG 
GGAAGGCAAAGGAAGATGGTG 


57 
57 


(ATC) 9 


FAM 


113-131 


1 


KF318156 


Lpi09 


F: 
R: 


AAATTTCCTCCGGCTGTTGG 
TCACTCGACGAATTTCCACC 


57 
57 


(ATC) 18 


PET 


254-316 


3 


KF318157 


LpilO 


F: 
R: 


AAGAGAGGTATGGGCGGAAC 
ACAGATTCGGAGTGGGAAGG 


57 
57 


(AAG) 9 


VIC 


236-283 


1 


KF318158 


Lpill 


F: 
R: 


CGTAATCTTCTGGCCTGCTG 
CCAGCTCCGGTATGATGTTG 


57 
57 


(CT) 9 


PET 


142-179 


1 


KF318159 


Lpil2 


F: 
R: 


GGGTGCTTTCGTTCCATTCC 
TGGTTTCATGGTGGAGAGGG 


57 
57 


(CT) n 


NED 


154-166 


1 


KF318160 


Lpil3 


F: 
R: 


ACAAAGGCCAGACAACAACC 
GCTGTGACTGTTGCTGTGAC 


57 
57 


(AGC) 9 


NED 


222-242 


2 


KF318161 


Lpil4 


F: 
R: 


CTTCCCAGGCAGTATCCCTC 
ATCGTGCTCTGTTCTACCGG 


57 
57 


(CTDn 


NED 


227-291 


3 


KF318162 


Lpil5 


F: 
R: 


GTTTGATAGAGCAGGCGACG 
CTATCGCGAATGACTGGCTTC 


57 
57 


(CTT) 9 


FAM 


98-119 


2 


KF318163 


Lpil6 


F: 
R: 


GAATCTTCTGCCCGCACTTC 
GGTAGCACTTTGAAGGCGTC 


57 
57 


(AGG) 12 


VIC 


158-194 


3 


KF318164 


Lpil7 


F: 
R: 


AATAGTCTCAGCCACTCCCG 
CGTCAAGCAGTTCGAATGGG 


57 
57 


(AT) n 


FAM 


309-351 


1 


KF318165 


Lpil8 


F: 
R: 


CCATGTTAAGGCACGGGAAC 
CCTACTGAGCCGGTGTACTC 


57 
57 


(ACG) n 


FAM 


374-380 


3 


KF318166 



Note: T a = annealing temperature. 

a Multiplex indicates loci that were mixed in the same capillary electrophoresis run. 
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